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Abstract. We present results from our magnetohydrodynamical simulations of ac- 
cretion flows onto black holes. Our main focus is the interplay between inflows and 
related outflows. We consider applications of such flows to the Galatic center and low 
luminosity active galactic nuclei. 



1 Introduction 

Dynamical evidence suggests that the nonhiminous matter within 0.015 pc of 
the Galactic center has a mass of sa 2.6 x 10 6 M Q (e.g., This matter 

is associated with Sgr A*, a bright, compact radio source [3], and provides very 
compelling evidence for the existence of a supermassive black hole (SMBH). 
Observations of Sgr A* in X-ray and radio bands reveal a luminosity substantially 
below the Eddington limit, ^Edd = 3 x 10 44 erg s _1 . For example, Chandra 
observations show a luminosity in 2-10 keV X-rays of ss 2 x 10 33 erg s" 1 (ten 
orders of magnitude below LEdd) Q • Chandra observations also revealed an X- 
ray flare rapidly rising to a level about 45 times as large, lasting for only ~ 10 4 s, 
indicating that the flare must originate near the black hole PQ . 

It is thought that the Sgr A* radiation is due to gas accretion onto the 
SMBH. Estimates for the accretion luminosity, L, rely on assumptions about 
the mass accretion rate, M a , and the efficiency of transforming the gas energy 
into radiation, n (i.e., L = nc 2 M a ). Both M a and 77 are uncertain and there is 
no generally accepted model which could explain the low luminosity of Sgr A* 
by predicting low enough M a or 77, or both. For example, Coker & Melia |S] 
estimated a rate of 10 -4 M© yr _1 from Bondi-Hoyle accretion of winds from 
nearby stars, while Quataert, Narayan, & Reid ^7] estimated the Bondi capture 
rate of 3 x 10 -5 M© yr _1 . On the other hand, the best fit spectral models of ^Sj 
and 19 have M a = 10~ 10 M Q yr" 1 and 10 -8 M Q yr -1 , respectively. Quataert, 
Narayan, & Reid J7| argue that the low luminosity requires M a to substantially 
sub-Eddington and sub-Bondi at large radius. However, the rate at which mass 
is captured into the accretion flow at large radii, the mass supply rate, does not 
necessarily have to be the same as the rate at which mass is accreted onto a 
black hole. 

Many aspects of spectral models for Sgr A* follow from the simple scaling 
laws of black hole accretion and will be present in any model, regardless of the 
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Fig. 1. Late time evolution of the mass accretion rate in units of the Bondi rate, for 
the fiducial model of MHD accretion flow in [15] (their model D). Time is in units 
of the Keplerian orbital time at the Bondi radius, which for this simulation was set 
at 1000 times the black hole radius. The figure shows in detail M a as a function of 
time toward the end of the simulation. Vertical arrows mark times corresponding to 
three (out of four) generic states of accretion: accretion is dominated by low-Z material 
which managed to reach the inner boundary despite a blocking corona and outflow 
from the torus (arrow d), accretion occurs only through the torus (arrow b) and no 
torus accretion but only very weak accretion through a very low density magnetized 
polar cylinder (arrow c). See Figs. 2 and 3 for the density maps and the velocity fields 
of the inner flow corresponding to the times marked. 



detailed dynamics. For a steady state model, an impressive range of spectra 
can be generated with the adjustment of several free parameters: the ratio of 
electron to ion temperature, the magnetic pressure, the radial density profile, 
the mass-accretion rate, and the mass-loss rate (e.g., JH])- Thus it is difficult 
to tightly constrain parameters of a steady state model. However, our recent 
magnetohydrodynamical (MHD) simulations present a way of constraining the 
accretion flow models by studying time variability [15] . As we will describe in 
next section, the simulations show that the inner accretion flow quasi-periodically 
changes both quantitatively (e.g., M a changes by ^ 1 order of magnitude) and 
qualitatively (e.g., the inner flow can be dominated by an equatorial accretion 
torus or a polar inflow). 



2 Results From MHD Simulations Of Accretion Flows 

We have performed axisymmetric, MHD, time-dependent simulations of slightly 
rotating accretion flows onto black holes [15] . Our simulations are complemen- 
tary to previous MHD simulations which considered strongly rotating accretion 
flows and started from a rotationally supported torus (e.g., [5], 0). We 



Accretion of Low Angular Momentum Material onto Black Holes 3 

attempt to mimic the boundary conditions of classic Bondi accretion flows as 
modified by the introduction of a small, latitude-dependent angular momentum 
at the outer boundary, a pseudo-Newtonian gravitational potential and weak 
poloidal magnetic fields. A weak radial magnetic field and distribution of the 
specific angular momentum, I, with latitude allow the density distribution at in- 
finity to approach spherical symmetry. Therefore our outer boundary conditions 
are consistent with X-ray images taken with Chandra which show that the gas 
distribution in the vicinity of SMBHs is close to spherical. 

We find that the material with high I forms an equatorial torus which can 
accrete onto the black hole because of magnetorotational instability (MRI) . The 
magnetized torus produces a corona and an outflow. The latter two can be 
strong enough to prevent accretion of low-? material through the polar regions 
(the source of accretion in the hydrodynamical inviscid case |14|h We find that 
the net mass accretion rate through the torus is lower than the Bondi rate and 
also lower than M a in the HD inviscid case. 

The accreting torus is the crucial component of our accretion flow. After 
an initial transient, the inner flow usually consists of a turbulent, gas pressure- 
dominated MHD torus with an outflow or corona or both, bounded by a magnetic 
pressure-dominated polar flow. The accretion through the torus is highly vari- 
able. In fact, it could be stopped for a relatively short time, by the accumulated 
poloidal magnetic flux that builds up during accretion and truncates the torus. 
However, we observe that even when the torus is truncated, there is inflow of 
material inside the torus and its mass and pressure build up. Consequently, the 
magnetic field is quickly pushed inward by the torus and the gas from the torus 
can again fall onto the black hole (note short-lived 'dips' and 'spikes' in the time 
evolution of M a shown in Fig. 1). Rapid time variability appears to be typical 
for MHD turbulent torii. Generally, we find that the properties of the inner flow 
(e.g., the radial density profile in the torus, properties of the corona and outflow) 
are very similar to those presented by |18j and |10| despite using different initial 
and outer boundary conditions. 

What is new in our simulations is the fact that the torus accretion can be 
supplemented or even replaced by stream-like accretion of the low-? material 
occurring outside the torus (e.g., Figs. 1, 2 and 3). When this happens, M a 
sharply increases and gradually decreases in a quasi-periodic manner. M a due to 
this 'off torus' accretion can be one order of magnitude larger than that due to the 
torus (Fig. 1). The off torus accretion is a consequence of our outer boundary 
and initial conditions which introduce the low-/ material to the system. This 
material can reach the black hole because the torus corona and outflow are not 
always strong enough to push it away. 

One can expect that in the vicinity of SMBHs some gas has very little angular 
momentum and could be directly accreted. Such a situation occurs likely at the 
Galactic center where a cluster of young, massive stars losing mass surrounds a 
SMBH. We propose then that the variability observed in Sgr A* (i.e., flares) could 
be due to the interplay between gas with a range of specific angular momentum 
similar to that found in our simulations. 
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3 Conclusions 

Our main results can be summarized as follows: 

• The properties of the accretion flow depend on an equatorial torus. 

• Accretion can be via the torus due to MRI and via the polar funnel where 
material has zero or low angular momentum. 

• The torus outflow and corona are natural mechanisms to narrow or even 
totally close the polar funnel for the accretion of low-^ material. 

• The net accretion rate for the MHD, is lower than for the Bondi flow and 
the HD inviscid flow. 

• Time variability of the inner flow may explain light curves in Sgr A*. In 
particular, we propose that the X-ray flares can be explained by quasi-periodic, 
off-torus accretion of the low I material. 

We acknowledges support from NASA LTSA grant NAG5-11736. We also 
acknowledge support provided by NASA through grant HST-AR-09947 from 
the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Fig. 2. Maps of logarithmic density overplotted by the direction of the poloidal velocity. 
This figure compares the inner flow in four different accretion states in our MHD 
simulations (Fig. 8 in 1151 1. We express the length scale in units of the black hole 
radius. The top left panel presents the two-dimensional structure near the beginning 
of simulations at t = 0.22 (not shown in Fig. 1). At this time, accretion onto the black 
hole occurs through both the torus and the polar funnel. The top right panel presents 
an example of an inner flow where accretion occurs only through the torus at t = 2.39 
(marked by arrow b in Fig. 1). The bottom left panel presents an example of an inner 
flow where there is no torus accretion but only very weak accretion through a very low 
density magnetized polar cylinder at t = 2.41 (marked by arrow c in Fig. 1, note that 
this state is very short-lived). Finally, the bottom right panel presents an example of 
an inner flow where accretion is dominated by low-Z material (i.e., the inward stream 
of gas directed present in the lower half of the panel). This material managed to reach 
the inner boundary despite a blocking corona and outflow from the torus at t — 2.36 
(marked by arrow d in Fig. 1). 
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Fig. 3. As Fig. 2 but with the r' and z' ranges increased by a factor of 5. Note the 
large scale, polar double-sided outflow from the thick equatorial torus in the top right 
and bottom left panels. The bottom right panel shows an accretion state where the 
outflow below the equator is replaced by an inflow. We propose that switching between 
this state and that presented in the top right panel (i.e., in Fig. 1, the states marked 
by arrows d and b, respectively) can explain radiation flaring in Sgr A*. 



